Electrons in an annealed environment: 
A special case of the interacting electron problem 
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The problem of noninteracting electrons in the presence of annealed magnetic disorder, in addition 
to nonmagnetic quenched disorder, is considered. It is shown that the proper physical interpretation 
of this model is one of electrons interacting via a potential that is long-ranged in time, and that its 
technical analysis by means of renormalization group techniques must also be done in analogy to 
the interacting problem. As a result, and contrary to previous claims, the model does not simply 
describe a metal-insulator transition in d — 2+e (e <C 1) dimensions. Rather, it describes a transition 
to a ferromagnetic state that, as a function of the disorder, precedes the metal-insulator transition 
close to d = 2. In d — 3, a transition from a paramagnetic metal to a paramagnetic insulator is 
possible. 



PACS numbers: 75.10.Lp; 71.30.+h; 64.60. Ak 
I. INTRODUCTION 

Local magnetic moments are known to play an impor- 
tant role in the behavior of disordered electronic systems, 
but the precise nature of that role remains incompletely 
understood. 1 One way to think about such local moments 
is that, in a disordered environment, the exchange inter- 
action between the electrons may be locally enhanced to 
the point where the electron spins order magnetically in 
a finite region in space. 2 The resulting magnetized re- 
gions are often referred to as local moments, or droplets, 
or rare regions. Since they are self-generated by the elec- 
tron system, they are in thermodynamic equilibrium with 
the other electronic degrees of freedom. It is therefore in- 
tuitively plausible that such local moments can be mod- 
eled as annealed magnetic disorder, in addition to the 
underlying quenched disorder that produces them. In 
Ref. 4 an explicit derivation has been given that corrob- 
orates this argument. There is experimental evidence 
for such local moments to influence the transport prop- 
erties of the electron system in important ways, and in 
particular they are suspected to influence the critical be- 
havior near the metal-insulator transition (MIT) that is 
observed in disordered electron systems. 1 However, the- 
oretically understanding the coupling between local mo- 
ments and transport properties has proven to be very 
hard. Studying and understanding the annealed disorder 
model mentioned above is expected to shed light on this 
important problem. 

Reference 4 provided such an analysis, and concluded 
that the annealed disorder leads to a new and very 
interesting type of MIT. The most exciting feature 
was that the transition was driven by the vanishing of 
the thermodynamic density susceptibility 9n/<9/i, and 
thus resembled a Mott transition more than an Ander- 



son transition. 5 This was even more surprising as the 
Coulomb interaction between the electrons, which is what 
usually causes a Mott transition, had not been explicitly 
taken into account in the model. 

Subsequently, Ref. 6 developed a general classification 
of quantum phase transitions with respect to, (1) whether 
one can describe the transition by means of a local order 
parameter, and (2) whether the order parameter suscep- 
tibility in the disordered phase is an analytic function 
of the wavenumber. The second criterion has an impor- 
tant bearing on which observables can become critical at 
a MIT: Criticality in d > 2 (d = 2 is the lower critical 
dimension for all known MITs of disordered interacting 
electrons) implies a logarithmic dependence on the renor- 
malization group (RG) length rescaling factor, and hence 
on the wavenumber, in perturbation theory in d = 2. 
This in turn implies a (weaker) nonanalytic wavenumber 
dependence in d > 2 away from criticality. 7,8 Although 
the considerations in Ref. 6 do not provide a mathemat- 
ically rigorous proof, they strongly suggest that dn/d[i 
cannot be critical at a MIT for a large class of models, 
which includes the model studied in Ref. 4. 

In the current paper we provide a thorough re-analysis 
of the model derived and motivated in Ref. 4, and re- 
solve this contradiction. We show that the RG analysis 
of the model performed in Ref. 4 had an incorrect struc- 
ture and led to unreliable results. A proper analysis of 
the model's renormalizability, and the resulting RG flow 
equations, show that dn/d/i is not singularly renormal- 
ized and hence not critical, in agreement with Ref. 6. In 
addition, it reveals that within a controlled e-expansion 
about d = 2, the model does not simply describe a metal- 
insulator transition. Rather, it displays a variant of the 
phase transition sequence that is known to occur in a 
related model with both quenched disorder and electron- 
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electron interactions (but no annealed disorder). 10 That 
is, as the disorder increases, there is first a transition 
to a ferromagnetic metallic state, and then, with further 
increasing disorder, a transition to a ferromagnetic in- 
sulator state. For d = 3 a transition directly from a 
paramagnetic metal to a paramagnetic insulator is pos- 
sible. 

This paper is organized as follows. In the next sec- 
tion, we give intuitive physical arguments that explain 
our model and our procedure to analyze it, and we sum- 
marize our results. In Sec. Ill we formally define the 
model and write it in a way that facilitates a renormal- 
ization group analysis. Sec. IV performs the renormaliza- 
tion to one-loop order, and Sec. V analyzes the results. 
Some technical issues regarding the model's renormaliza- 
tion properties are relegated to Appendices A and B, the 
flow equations for the interacting and annealed disorder 
models are compared in Appendix C, and a perturbative 
analysis of the free energy is given in Appendix D. 

II. PHYSICAL ARGUMENTS 

Since some of our detailed arguments are quite techni- 
cal, we start by giving some intuitive physical arguments 
to explain both our general strategy and our results. 

A. Annealed disorder as a model for local moments 

We start by recalling the argument for why annealed 
disorder models local moments. 4 Any field theoreti- 
cal treatment of a statistical mechanics problem starts 
with a functional integral representation of the partition 
function, 1 

Z = J D[<j>] e- s[ * ] . (2.1a) 

The form of the action S defines the model under con- 
sideration, and the mathematical nature of the field <fi 
depends on whether the system is classical or quantum 
mechanical, consists of fermions or bosons, and whether 
the model is a microscopic one in terms of fundamental 
fields, or of an effective nature. The usual procedure is to 
identify a saddle point of S that approximately contains 
the physics one is interested in, to expand about this 
saddle point, and to employ perturbation theory and the 
renormalization group. In a system with quenched disor- 
der there will be, apart from homogeneous saddle-point 
solutions, solutions where the field <f>, or some compo- 
nents of it, are nonzero only in certain regions in space. 
Such inhomogeneous saddle points have been proposed 
as a description of rare regions in classical magnets by 
Dotsenko et al. 2 This concept was generalized to quan- 
tum magnets, 3 and to the effective field theories used to 
describe MITs in quenched disordered electron systems. 4 
In a large system there will be many rare regions that 



interact only very weakly, and thus exponentially many 
almost degenerate saddle points, since the orientation of 
the field on the rare regions is arbitrary. These saddle 
points are expected to be separated by large energy bar- 
riers, and thus to not be perturbatively accessible from 
one another. Within perturbation theory, and denoting 
the n-th saddle-point field configuration by and the 
fluctuations by ip, one can therefore write the partition 
function 

Z^fl[^- s|$W ^ . (2.1b) 

n 

In the thermodynamic limit, the discrete set of saddle 
points becomes a saddle-point manifold that needs to be 
integrated over. The saddle-point field configurations $ 
thus become degrees of freedom that are governed by 
some probability distribution P[<&], are integrated over 
at the level of the partition function, and couple to the 
field tp by means of some coupling S c that is determined 
by the action S, 

J £>[$] P[$] J D[<p] e -SM+Sc[*,y] _ ( 2 .lc) 

They therefore act like annealed disorder. Note that 
in giving Eq. (2.1c) we implicitly assume that the (an- 
nealed) disorder adjusts and comes to equilibrium with 
the fluctuations ip. If the disorder were fixed on the time 
scale of the ^-fluctuations, then it would be quenched 
disorder. In the latter case, for the average over saddle- 
points to be meaningful, \nZ rather than Z should be 
averaged over the <I>-fields. 12 

In our case we are interested in rare regions that carry 
a magnetic moment. According to the arguments recalled 
above, they can be modeled by annealed magnetic disor- 
der in addition to the quenched disorder that allows for 
the inhomogeneous saddle-point solutions. In the sim- 
plest possible model the annealed disorder has a Gaus- 
sian distribution, and is static. The latter means that 
the coupling constant, or the annealed magnetic disorder 
strength, will be proportional to the temperature. 4 This 
is just the Boltzmann weight assigned to these classical 
degrees of freedom that are in equilibrium with the elec- 
trons. We emphasize that this model, and its derivation 
in Ref. 4, is unaffected by our considerations concerning 
its analysis and interpretation, which differ from the one 
given in that reference. 

B. Annealed disorder as an effective interaction 

The physical effects of annealed disorder are funda- 
mentally different from those of quenched, or frozen-in, 
disorder. 12 The former gets integrated over at the level 
of the partition function, cf. Eq. (2.1c), the latter, at 
the level of the free energy. Consequently, integrating 
out annealed disorder generates a physical effective in- 
teraction between the degrees of freedom that couple to 
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it, the electron spin density in our case, which can be 
understood as resulting from an exchange of annealed 
disorder fluctuations between the electrons. The effects 
of quenched disorder, on the other hand, are more subtle 
and fundamentally different from those of interactions. 

It is therefore plausible that a system of noninteract- 
ing electrons in the presence of both quenched and an- 
nealed disorder will behave in many respects like one 
with quenched disorder only and an additional electron- 
electron interaction. As the only difference one would 
expect that, if the annealed disorder is modeled as static, 
the resulting effective interaction will be infinitely long- 
ranged in time, a feature that one would not expect to 
have qualitative effects. This expectation is in contra- 
dition with the results of Ref. 4, which found behavior 
that was drastically different from that of electrons in- 
teracting via an instantaneous interaction. In particular, 
this reference predicted a MIT of Mott type, where the 
thermodynamic susceptibility dn/d/j, vanishes. This is in 
contradiction to both explicit calculations for quenched 
disordered, interacting electron systems, which find that 
dn/d/j, is not singularly renormalized, 13,1 and very gen- 
eral considerations in Ref. 6. 

The analysis that will be presented below removes this 
contradiction, and illustrates the technical issues behind 
the above intuitive physical considerations. We will show 
that the technical treatment of the annealed disorder in 
analogy to that of quenched disorder in Ref. 4 was not 
only in disagreement with the above physical arguments, 
but led to an unnatural structure of the theory. This 
in turn led to incorrect assumptions about the behavior 
under renormalization, and ultimately to physically in- 
correct results. A treatment of the annealed disorder in 
analogy to an interaction, on the other hand, does not 
run into these problems and yields results that are in 
agreement with all known constraints. 

III. THE MODEL AND ITS 
RENORMALIZ ABILITY 

In this section we consider the same effective field the- 
ory as in Ref. 4. 

A. Effective field theory 

Our starting point, as in Ref. 4, is Wegner's nonlinear 
sigma- model (NLcrM) 14 for noninteracting electrons with 
nonmagnetic quenched disorder. The action reads 

-Anl.m - 7^J dxtT [VQ(x)] 2 + 2tf«y"dxtr [flQ(x)]. 

(3-1) 

Here Q(x) is a matrix field that comprises two fermionic 
degrees of freedom. Accordingly, Q carries two fermionic 



Matsubara frequency indices n and m, and two replica 
indices a and (3 to deal with the quenched disorder. The 
matrix elements Q"m are spin-quaternion valued to allow 
for particle- hole and spin degrees of freedom. It is con- 
venient to expand them in a basis t t (g> Sj (r, i = 0, 1, 3.3) 
where to = sq is the 2x2 unit matrix, and 71,3.3 = 
— si,3.3 = —i<?i,3.3, with Oj the Pauli matrices, 15 

Qn& = £E&i • (3.2a) 

r i 

For simplicity, we will ignore the particle-particle or 
Cooper channel, which amounts to dropping n and r 2 
from the spin-quaternion basis. 15,1 The Q"m are then el- 
ements of C x Q, with C and Q the complex number field 
and the quaternion field, respectively. The l r Q nm obey 
the following symmetry properties (for r = 0, 3), 16 

rQnm = ( — ) rQmn ' (3.2b) 
i rQ a ni = (-Y +li rQ Z , (* = 1, 3.3) . (3.2c) 

Alternatively, we can write the spin indices explicitly, 
and consider matrix elements Qnm.ij tnat are complex 
number valued. Q is subject to the constraints 

Q 2 (x) = l , trQ(x) = . (3. 2d) 

These constraints are conveniently implemented by 
parametrizing Q in terms of matrices q whose matrix 
elements, g"^, are restricted to frequency labels n > 0, 
m < 0. In terms of the q, Q can be written in block 
matrix form 

V ? T -y/l-qlq+lj 

Here the block matrices, clockwise from the upper left, 
correspond to frequency labels n, m > 0; n > 0, m < 0; 
n, m < 0; and n < 0, m > 0, respectively. 

= 5 nm 5 ati Vl n (r <g> s ) in Eq. (3.1) is a frequency 
matrix with fi n = 27rTn a bosonic Matsubara frequency 
and T the temperature. G is a measure of the disorder 
that is proportional to the bare resistivity, and the fre- 
quency coupling ijW is proportional to the bare density 
of states at the Fermi level, tr denotes a trace over all 
discrete degrees of freedom that are not shown explicitly. 

The properties of this model are well known. 14 ' 15 ' 1 
The bare action describes diffusive electrons, with D = 
1/GH^> the diffusion coefficient. Under renormaliza- 
tion, D decreases with increasing disorder until a MIT is 
reached at a critical disorder value. The critical behav- 
ior is known in an e-expansion about the lower critical 
dimension d = 2. In the absence of the Cooper chan- 
nel, the MIT appears only at two-loop order at a criti- 
cal disorder strength of 0{y/e). which determines 
the specific heat coefficient, the spin susceptibility, and 
dn/dfi, is uncritical, which makes this MIT an Anderson 
transition. 
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Now we add magnetic annealed disorder to the model. 
The motivation for this is the fact that annealed disorder 
models certain types of local moments, see Sees. I and II 
above. A technical derivation of this has been given in 
Rcf. 4, and the main idea has been recapitulated in Sec. 
II A. Annealed disorder implies that the Q in the result- 
ing terms all carry the same replica index; 12 otherwise, 
the functional form of the resulting additional term in 
the action can be taken from Ref. 15, which considered 
quenched magnetic disorder. From that reference, we 
have 

A£1 = 2TJW /dx^^tr [(T 3 ® S ,)g Q «(x)] 2 . 

J a j=l 

(3.3a) 

The coupling constant JW is a measure of the strength 
of the magnetic disorder. The temperature prefactor in 
Eq. (3.3a) is a consequence of the static nature of the 
local moments considered within this model, as has been 
explained in Rcf. 4 and Sec. II above. Equation (3.3a) is 
the only annealed magnetic disorder term if fluctuations 
of the matrix field Q on all length scales are taken into 
account in calculating the partition function. However, 
the NLctM is an effective theory for long- wavelength fluc- 
tuations, and it is therefore convenient to project the an- 
nealed disorder term onto this regime as well. It has been 
discussed in detail in Ref. 16 that this can be achieved 
by means of a phase space decomposition and a relabel- 
ing of momenta. Applied to Eq. (3.3a), this procedure 
generates another contribution to the action, 




(3.3b) 

The coupling constant is in general independent of 

^ • A&nn 

and A^Sn enter the action additively with 
the understanding that only long-wavelength fluctuations 
are integrated over in calculating the partition function. 
Note that in the case of quenched magnetic disorder, a 
complete phase space decomposition leads to a term anal- 
ogous to Eq. (3.3b), but it is zero in the replica limit 
because the replica sum is then part of the trace. 

As we will see, under renormalization the annealed dis- 
order terms generate another contribution to the action 
that takes the form 

A^ = 2HW J rfx tr [sgntt Q(x)] , (3.4) 

so we add this right away. For a discussion on why this 
term must be present on physical grounds, see Section 
VB. 

A^A NLr7M +A^ n +A^ n +A^ , (3-5) 



is the complete action for our model, and the partition 
function is obtained as the functional integral 

Z = [ D[Q] S[Q 2 - 1] e A ^ . (3.6) 



B. Annealed disorder as a long-ranged interaction 

A = -4nl<tM + -4ann defines the model studied in Ref. 

(2) 

4. ^4ann was neglected in that reference, but this term 
will not be of crucial importance in what follows. Terms 
that appear under renormalization and indicate the ap- 

(2) 

pearance of A n were interpreted differently in Ref. 4, 
and we will discuss this point in Sec. VB below. A re- 
lated point is that we have written „4inn in a form that 
is different from the one in Ref. 4. The latter representa- 
tion was modeled after the way one would treat quenched 
disorder, and it added and subtracted a term where all 
replica indices of the Q are not the same. As we will see, 
this formulation, which is a matter of taste at this point, 
is rather unnatural at the stage of a RG analysis, and 
this led to the incorrect RG treatment of the model in 
Rcf. 4. We therefore write the annealed disorder term in 
a form that is strictly diagonal in the replica index. This 
replica structure is common to both the annealed disor- 
der term, and any electron-electron interaction term, and 
one would therefore expect the renormalization proper- 
ties of the current model and one of interacting electrons 
to have common features. To underscore this point, we 
rewrite the annealed disorder part of the action by split- 
ting it into spin-singlet and spin-triplet contributions, 

with 

4 1 nn ) = Z ^ (1 ' S) £££(-r 

nm a r— 0,3 

xtr [(r r ® s ) (x)] tr [(r r ® s ) Q« a „(x)] , 

(3.7b) 

^ = jj (1,t) EEE(-) r E 

nm a. r— 0,3 i—1 

xtr [(r r ® 8i ) Q™ (x)] tr [(tv ® 8i ) Qrjx)] , 

(3.7c) 

>e> = EE E(-) r E 

nm a. r— 0,3 i—1 
xtr [(TV ® 8i ) QZ(x)\ tr [(TV ® Si) QZnW] , 

(3.7d) 

where we have used Eqs. (3.2a) - (3.2c). Here 
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and 



J( 2,t) = 8J (2) /7r 



(3.7e) 



(3.7f) 



This relation between the bare values of J^ 1 '^ and j( 1,t; ) 

will be important later. Notice that J^ 1 ^ < 0, while 
j(i,t) > Q) j(2,t) > 

Comparing these expression to the correspoding ones 
for an electron-electron interaction, 16 one sees that they 
have the same structure except for the frequency sector. 
Transforming from Matsubara frequency space into time 
space reveals that the annealed disorder corresponds to 
an interaction that is infinitely long-ranged in time. This 
is physically plausible, as has been explained in Sec. II B. 



C. Renormalizability considerations 

For reasons explained in Appendices A and B, we will 
choose a field-theoretic RG method 9 over a momentum- 
shell RG. 17 Before we start analyzing our model by means 
of this method, we need to ask whether the model is 
rcnormalizable, and how many renormalization constants 
are required. Much is known about the renormalization 
properties of the NLctM, Eq. (3.1), with additional in- 
stantaneous interaction terms. The pure NLctM is known 
to be renormalizable with two renormalization constants, 
one for the coupling constant G and one field renormal- 
ization constant. 9 The frequency coupling ijW turns out 
to not carry a renormalization constant of its own. In 
the presence of an instantaneous interaction, the proof 
of renormalizability for the NLerM breaks down, and the 
renormalizability of the model has never been proven. 
However, there is much evidence that the model is still 
renormalizable, with two additional renormalization con- 
stants for the interaction, and with acquiring a 
renormalization constant of its own. The two renormal- 
ization constants for the interaction terms correspond to 
symmetric and antisymmetric combinations of terms bi- 
linear in Q, respectively. 18 ' 19 The same arguments ap- 
ply to the present model, and are given in Appendix 
B. From Eqs. (Bl), we conclude that we need to write 



^ann 



A™ 



A^l ' s \ and analogously split A^nn and 



-4inn\ With 



J r,m n n (1 Q 



a a 
nm 



(x) 



nm a r— 0,3 



A^ = -2,TJ^j dx££ £ [?Ox) 



(3.8a) 



nm a r— 0,3 



(3.8b) 



A^ = -2nTJ** /dxEE E [E^:» 



nm a r— 0,3 i—1 

3 \ 

— 1 i /-i QQ / \ i r\ aa I \ 



(3.8c) 



A^ = -2nTJ^ /dxX)E E [Er«W 

" . — „. n o „• i 



nm a r— 0,3 i—1 
3 




X rC( X )-oE( - ) rCnWfCW 



(3.8d) 



nm a r— 0,3 i—1 



(3.8e) 



A™ = 2 7 rTji 2 ' t) 



/ dx EE E [E : 

nm a r— 0,3 i—1 
/ 3 \ 



<5nn(x) 



X^r m (x)-^E( I ) rQZ^VrQZ^) 



(3.8f) 

In writing Eqs. (3.8) we have made use of Eqs. (3.2a) 



(3.2c). The symbol 



is a shorthand for 35m — 



Y^j=i&ij- The J^'^ are coupling constants whose bare 
values are equal, 

J^' s) = Ji 1,s) = J^ s ) , (3.9a) 
but in general they renormalizc differently. Similarly, 

= j(i,t) = J( i,t) ; (39b) 

(3.9c) 



j(2,t) = j(2,t) = j(2,t) 



in the bare theory, but under renormalization these 
equalities do not in general remain valid. All of the J + 
require only one renormalization constant, which we will 
denote by Z + , and the J_ require another one, Z-. In 
addition, a renormalization constant for is needed. 

In addition to the relations given by Eqs. (3.9), there is 
the relation between J^ 1 ^' and given by Eq. (3.7f). 

It will turn out that these constraints leads to a degen- 
eracy in the RG flow. This is most easily handled by 
relaxing the condition, Eq. (3.9a). Instead of Eqs. (3.9) 
and (3.7f) we therefore write 



j(l,s) = J(1 , B) ± A 



(3.10a) 
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j£,t) = j(i,t) ) (3 10b) 

j£ 2 ' t} = J( 2 -*) , (3.10c) 

and 

j(i.s) +3j(M) = o . (3.10d) 

Choosing A ^ will remove the degeneracy in the RG 
flow. In the end, we will consider the limit A — > to 
obtain the behavior of our original model. 



APf(k)=Cf(k)-P n (k) , (4.2f) 

and 

J s = \ (J| M) + J1 M) ) - \ (Ji 2,t) - , (4.2g) 

J t = l(j| 1 ' t) + ji 1 ' t) ) + ±(j< 2 ' t >-J< 2 ' t >) . (4.2h) 

2. One-loop corrections 



IV. RENORMALIZATION TO ONE-LOOP 
ORDER 

A. Perturbation theory 

i. Gaussian propagators 

We now perform a one-loop RG analysis of the model 
defined in Sec. III. To this end, we expand the action in 
powers of the matrix q defined by Eq. (3.2e). To Gaussian 
order we find 

^^EEE r«i 2 (p) r(2) (p. n ni _„ a ) j.«i 2 (-p) . 



p 12 i,r 



(4.1a) 



with 



r< 2 )(k, n„) = k 2 /G + H^n n + r^H-kT 

+8 aia2 irTG [S i0 J s + (1 - S i0 )J t ] . (4.1b) 

the bare two-point vertex. The Gaussian g-propagators 
are obtained by inverting this quadratic form. We find 

G 

(r<?i 2 ( k )i9 3 4(p)> =^- P S 13 5 24 S rs S i: j - l V 12 (k) . 

(4.2a) 

Here (. . .) denotes a Gaussian average, and 1 = (ni, ai), 
etc., are indices that comprise both the Matsubara fre- 
quency index and the replica label. The propagators l T> 
read 

(k)+6 aia2 AV s ni _ n2 (k) , (4.2b) 
1 ' 3 - 3 P 12 (k) = V ni _„ 2 (k) + S aia2 AV' ni _ n2 (k) , (4.2c) 

where 

1 



^n(k) 
2?f(k) 



k 2 + GHWfl n + GH( 2 )2ttT ' ( 4 ' 2d ^ 
1 

k 2 + GHWn n + GH^2ttT + GJ^ttT ' 

(4.2e) 



By expanding the action to 0(q 4 ) and calculating all 
diagrams with the topological structure shown in Fig. 
1, we obtain the one-loop corrections SG, 5H^ X \ etc. to 




FIG. 1. Structure of diagrams that renormalize the 
two-point vertex. 

the coupling constants in the Gaussian propagators, Eqs. 
(4.2), or the 2-point vertex, Eq. (4.1b). The explicit cal- 
culation is similar to the one for the case of an instanta- 
neous interaction, 1 but substantially simpler due to the 
absence of cubic terms in the g-expansion. We hnd 



5G=^ (K + + K_)I 2 , 
5H W = Z^1 {K+ + K _ )l2 ; 



(4.3a) 



16 



(4.3b) 



-G 
16 



H^{K + +K_)+' A - (j™ + J™ 



J + 



(2,t) 



1 



J 



(2,t) 



2 j| 2 ^ + 2 



G 



SJ,= 



32 
-G 

1 



^U 2 + 3J 2 )7 2 + ^(jr + J 



r(2,t) 



-G 



L- I 



(M) 



(4.3c) 



(4.3d) 



+ 2 J + 2 J 



(2,t) 



G 



16 



j(2,t) + j(2.t) 



h 



with 



(4.3e) 
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Here we have defined linear combinations of coupling con- 
stants, 

K± = j£ 1>s) + 3 j£ M) = ±A , (4.4a) 
L± = j£ 1,s) - j£ M) , (4.4b) 

where the second equality in Eq. (4.4a) is due to Eqs. 
(3.10). This will be important later. We have also defined 
one-loop integrals 

h = G J dp V n {p) = -G/Ge , (4.5a) 
h = G [ dp 2irTj2(Vn(p)f = -G/GH^e . (4.5b) 

J n 

Here e = d — 2, and G = GSd/{2ir) d with Sd the surface 
area of the (d — l)-sphcre. In giving the second equalities 
in Eqs. (4.5) we have chosen to use dimensional regular- 
ization, and in what follows we will use a field-theoretic 
RG method. At a perturbative level, this is a matter of 
choice, and we could just as well use the momentum-shell 
RG method. In that case, the factors of — 1/e in Eqs. 
(4.5) would be replaced by In 6, with b the RG length 
rescaling factor. For arguments that go beyond pertur- 
bation theory, however, it is advantageous to use the field 
theory approach, as is explained in the Appendices A and 
B. 

In addition to these renormalizations of the two-point 
propagator or vertex function, we will also need the one- 
point vertex to one-loop order. This is given by the 
diagram shown in Fig. 2, and a simple calculation yields 




FIG. 2. Structure of diagrams that renormalize the 
one-point vertex. 



r (1) ^(K:(x))- 1 = i--(^ + + X-)/ 2 • (4.6) 

For later reference, we notice that the one-loop correc- 
tions to G, and vanish in the limit A — > 0, and 
that 

<5J s + 3<5J t =0 , (4.7) 

as can be seen by using Eqs. (3.10). Furthermore, a cal- 
culation shows that 

5H& +SJ s = ^Ah , (4.8) 

which also vanishes as A — > 0. 



B. Renormalization 

1. Renormalization constants 

We now proceed to renormalize the theory, i.e., we ab- 
sorb the singularities in the e — > limit that are present in 
perturbation theory into renormalization constants. We 
define renormalized coupling constants g, ft/ 1 ), etc., by 

G = vr*Z g g , = zWhW , = Z™h™ , 

KM _ 7 ,-(l,s) T(l,t) _ 7 -(l,t) j(2,t) _ 7 .(2,t) 
j(l,-) = Z _fi,*) ; j(l,t) = z _£l,t) ^ j(2.t) = z _.(2,t) j 

(4.9) 

where \x is an arbitrary momentum scale. The renormal- 
ization statement is 9 

r^p.fi,,; g, h,j+,j-;fj.) = 

zW 2 )rW(p,0„;G,^,J + ,J_) . (4.10) 

Here is the renormalized TV-point vertex function, 
Z is the field renormalization constant, and H and J± 
represent the various frequency and annealed disorder 
coupling constants. The assertion that all vertex func- 
tions can be made finite to all orders in the loop expan- 
sion by the five renormalization constants defined in Eq. 
(4.2), plus the field renormalization constant, is equiv- 
alent to saying that the theory is renormalizable with 
these renormalization constants. As we have mentioned 
before, there is strong evidence for this statement to be 
true, which is recapitulated in Appendix B, but it has 
not been rigorously proven. 

Assuming that the theory is renormalizable, the six 
equations, Eqs. (4.3) and (4.6), suffice to determine the 
six renormalization constants to one-loop order. While 
it is possible to do so for arbitrary bare values of the 
coupling constants, the results simplify substantially if 
one uses Eqs. (3.10). Using minimal subtraction, 9 and 
taking the limit A — ► 0, we obtain 



z 


= l + 0(g 2 ) , 


(4.11a) 


z g 


= l + 0{g 2 ) , 


(4.11b) 


7 (1) 


= l + 0(g 2 ) , 


(4.11c) 




= l + £ K (g,h,j+,j-)/hW , 


(4.11d) 


z+ 


g 20 s ( 3 ,ft,j+,j_) 2 
e + 7 -(i.») + U{9 ' ' 


(4.11c) 


z_ 




(4.11f) 



Here k and Sjt are functions of the renormalized coupling 
constants that are given by SH ^ and £J Sit as functions 
of the bare ones, 
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k(G,H, J) 
B ,t(G,ff, J) 



-e8H^(G, H, J)/G 
-e5J s , t {G,H, J)/G 



An inspection shows that, in the limit A 
n(G,H,J) = -<p Stt (G,H,J) . 



0, 



(4.12a) 
(4.12b) 



(4.12c) 



Notice that Z + = Z_, at least to one-loop order. Since 
the bare values of the various J± are identical, this means 
that the renormalized values are also identical, and we 
can drop the distinction between the j+ and the j- . We 

will thus write j+' s ^ = j-' s ^ = j^' s \ etc. We note that 
this is a consequence of the relations expressed by Eqs. 
(3.10), and would not necessarily be true for more general 
models. 



2. Flow equations and their solutions 

We now are in a position to determine the RG flow 
equations for the coupling constants. Defining I = — In fd 
(or £ = In b in an alternative momentum-shell approach) , 
and using Eq. (4.12c), we obtain from Eqs. (4.9) and 
(4.11), 

dg_ 

di 
dhW 
di 

di 
dj^ 
dl 

The flow of the remaining coupling constants j can be 
obtained by relating them to j^ 1 '*'. This is a consequence 
of there being only two renormalization constants for all 
of the J. We obtain 





(4.13a) 


0(g 2 ) , 


(4.13b) 


-gcj) s (g,h,j) + 0(g 2 ) , 


(4.13c) 


^-<j, s (g,h,j) + 0(g 2 ) . 


(4.13d) 



J 



(M) 



;(2,t) 



■{i,t) j(i.s) /j(i,t) = _3 j(M) 
J(M) 



;(M) 



(4.13e) 
(4.13f) 



In order to determine the nature of these flows, we 
calculate 4> s from Eqs. (4.3d) and (4.12b). We find 



-3 (j (M) ) 2 

M g, h , j ) = —— W - 



1 



.7 



[2d 



4(j(i>*))2 



0(9) 



(4.14) 



We see that 4> s < 0, unless j( 2 '*) is larger than (J^ 1 '*)) 2 
in suitable units (note that the J's and iJ's all have the 
dimensions of a density of states). This makes physical 
sense: From Eqs. (3.7c) and (3.7d) we see that and 
_4(2,t) are S pi n _triplet interactions with different signs. 
j(M) > q p romo tes ferromagnetism, and J^ 2 '*' > 
weakens that tendency. In two-dimensions, for physi- 
cally sensible values of the coupling constants, we thus 



have 4> s < 0, and and j^ 1 ''-* both scale to infinity. In 
d > 2, the RG flow equations can be solved explicitly and 
shown to describe a quantum phase transition by intro- 
ducing, as in Rcf. 10, a scaling variable y = gj^'^/h^ 
that obeys 



f e= -ey + y y2 + 0(y 3 ) 



(4.15) 



We see that Eq. (4.15) allows for a fixed point value y* = 
e + 0(e 2 ). Denoting the deviation from this fixed point 
value by Sy, we find 



and 



5y(b) = 5y(b=l)b^°^ , 
h^(b) = h^(b=l)b^ +0 ^ , 

j(U) { p )=j (U) { p =l)b e+0(e*) 

h^(b) = h^(b=l)b° + °^ , 
g{b)=g{b=l)b-*+°^ . 



(4.16a) 

(4.16b) 
(4.16c) 
(4.16d) 
(4.16e) 



The behavior of all observables of interest can be deduced 
from the above flows, see Sec. VA below. 



V. DISCUSSION 

A. Physical interpretation, and results 

For a physical interpretation of our results we first need 
to relate physical observables to the coupling constants 
of our theory. Some observables can be identified directly 
in analogy to the corresponding identification in the case 
of an instantaneous electron-electron interaction. From 
the derivation of the NLctM, G is known to be related to 
the bare conductivity a via 14 ' 1 



/ttG . 



(5.1a) 



The single-particle or tunneling density of states N at an 
energy u from the Fermi level is related to the one-point 
vertex by 13 

N(e F +uj) = ^ (T (1) ) (iw n -» w + iO) . (5.1b) 

Equations (4.13a) and (4.11a) show that a and N are not 
renormalized, at least to one-loop order, 



da 

~d£ 
dN 

~d£ 



= 0( 5 2 ) , 
= 0( 3 2 ) . 



(5.2a) 
(5.2b) 



The scaling behavior of the relevant operator Sy, Eq. 
(4.16a), determines the correlation length exponent. De- 
noting the dimensionless distance from the critical point 
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by t, and the correlation length by £, one finds for small 
t 

e«|tr v , (5.3a) 

with a correlation length exponent 

v = l/e + 0(l) . (5.3b) 

Other quantities of interest are various susceptibilities, 
in particular the specific heat coefficient jy — Cy/T, 
the spin susceptibility Xs, and the density susceptibility 
dn/dfi. Their relations to the coupling constants in the 
field theory are less obvious. We therefore use scaling 
arguments, in conjunction with perturbation theory for 
the free energy, to determined their respective cricital 
behavior. We start with a homogeneity law for the free 
energy. From the Gaussian propagators, Eqs. (4.2), we 
see that that in principle there are three different time 
scales in the theory, given by 

ti = C d 9h {1) ~ T 1 , (5.4a) 

r 2 = ev (2) ~r 2 , (5.4b) 

T3 = i d 9j {U) ~ Z" , (5.4c) 

Here 21,2,3 are the dynamical exponents related to these 
time scales. To one-loop order we have 

21 = d-e + 0{e 2 ) = 2 + 0(e 2 ) , (5.5a) 

22 = 23 = d-e + e + 0{e 2 ) = d + 0(e 2 ) , (5.5b) 

leaving us with two times scales and dynamical expo- 
nents. The free energy density / therefore has two dif- 
ferent scaling parts, and we can write 

f(t, T, . . .) = b- {d+z ^ h{t b 1 ^, T b Zl ,Tb Z2 ,...) 

+b -(d+z2) f 2 (tb 1 l v ,Tb t \TV*,...) . (5.6) 

Here /1 and / 2 are scaling functions, and the ellipses 
denote the dependence of / on external fields that are 
not shown explicitly. 

The specific heat coefficient is obtained by differenting 
/ twice with respect to T. The leading contribution is 
obtained by differentiating f± with respect to the tem- 
perature scale that carries the dynamical exponent z 2 . 
This yields 

-Yv(t)~\t\- a , (5.7a) 

with a critical exponent 

a = v(2z 2 -d-z 1 ) = l + 0(e) . (5.7b) 

To ascertain that this leading contribution has a nonzero 
prefactor we check against perturbation theory for the 
free energy, which is given in Appendix D. From Eqs. 
(Dlb,Dlc) we see that there is indeed a contribution 
from differentiating twice with respect to the tempera- 
ture in the propagators, which carries a dynamical ex- 
ponent 22. The temperature prefactor in the expression 



/ = —{T/V)h\Z for the free energy density has been 
absorbed into the frequency integration measure. The 
frequency, however, scales like a wavenumbcr squared, 
and therefore carries an exponent 21 . 

A very similar argument applies to the spin suscepti- 
bility. A magnetic field B couples to the electrons via a 
Zeeman term (amongst other coupling mechanisms), and 
hence can scale like an energy or temperature. The spin 
susceptibility is obtained by differentiating / twice with 
respect to B, and once therefore expects Xs to scale like 
the specific heat coefficient, viz. 

Xs(t)~|*r 7 . (5-8a) 

with a critical exponent 

7 = a = l + 0(e) . (5.8b) 

Again, perturbation theory confirms that the leading 
contribution obtained in this way is nonzero. This is 
easily seen from Eqs. (Dlb,Dlc) by taking into account 
that B 7^ leads to a mass /zb-B in two of the spin-triplet 
propagators that contribute to the Gaussian approxima- 
tion for the free energy. 

Finally, we consider dn/d/j,. Although the chemical 
potential \i is dimensionally an energy, it differs funda- 
mentally from either T or /j-bB, since it represents the 
microscopic energy or inverse time scale. As such, it must 
have an effective scale dimension of zero. Consequently, 
we obtain from Eq. (5.6), by differentiating twice with 
respect to /j, 

(dn/dfj,)(t) = const. + 0(f<* fzi >) . (5.9) 

dn/dji thus has only a weak nonanalytic i-dependence in 
addition to a leading noncritical contribution. Again, 
this is consistent with perturbation theory: The only 
/i-dependence of the free energy, Eq. (Dlb), is through 
the various coupling constants in the propagators. All 
of these multiply either a frequency or a temperature. 
Differentiation with respect to /i therefore does not pro- 
duce a singular integral unless / itself becomes singular. 
Power counting shows that this happens only for dimen- 
sions d < —2, in agreement with Eq. (5.9). This failure of 
differentiation with respect to a field to produce a singu- 
larity is an illustration of a more general argument given 
in Ref. 6. 

The physical interpretation of these results is now 
clear. The RG flow at one-loop order is qualitatively the 
same as for electrons interacting via an instantaneous 
interaction, see the comparison between the two flows 
given in Appendix C. In the latter case, the runaway 
flow of the equivalent of j^ 1 '*-* (fc t in Appendix C) at 
one-loop order in d = 2 suggests a ferromagnetic ground 
state. In d = 2 + e there is a phase transition where 
the homogeneous magnetic susceptibility diverges. This 
transition has been identified with a ferromagnetic phase 
transition where the magnetic susceptibility diverges like 
Xs ~ |*| — K , as in Eq. (5.8a). 10 ' 20 The runaway flow thus 
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simply reflects the fact that t is RG relevant at a fer- 
romagnetic transition. The result of this interpretation 
agrees with a more direct, and more explicit, theory for 
the ferromagnetic transition. 20 In the current case, the 
theory describes an infinite-range version of this transi- 
tion, due to the interaction being infinitely long-ranged 
in time. These considerations strongly suggest that the 
physical results we have derived above to one-loop order 
actually hold to all orders in the loop expansion, as they 
do in the instantaneous interaction case. 10 In particular, 
we expect that dn/d^i is not rcnormalized to all orders, 
in agreement with Ref. 6. It also follows that the phase 
diagram for the present model is qualitatively similar to 
the one for the interacting case, with a ferromagnetic 
transition always preceding an MIT for d > 2, while for 
d = 3 a direct transition from a paramagnetic metal to 
a paramagnetic insulator is possible. 10 ' 1 There are, how- 
ever, differences in the detailed properties of the transi- 
tion as compared to the one studied in Refs. 20 and 10. 
For instance, in the latter the specific heat has a much 
weaker singularity than the spin susceptibility, while here 
they show the same scaling behavior. In this respect the 
current case is reminiscient of the Brinkman-Rice theory 
of the Hubbard MIT. 21 

Although the transition in the present model is clas- 
sical, in the sense that the order parameter is purely 
static, it couples to quantum mechanical degrees of free- 
dom in the form of the diffusive electrons. An explicit 
description of the transition could be obtained along the 
lines of Ref. 20. However, given the schematic nature 
of our model, we will not pursue this here. The same 
conclusion, namely that the model under consideration 
describes a ferromagnetic transition of a classical nature, 
has recently been reached by Vojta and Narayanan by 
means of very different arguments. 22 We stress, however, 
that our goal here has not been to describe a magnetic 
transition. Rather, it was to resolve the conflict between 
the results of Refs. 4 and 6, and to check whether or not 
our model of electrons with both quenched and annealed 
disorder describes an unusual MIT. As we have seen, the 
answer to the latter question is negative. 

B. Comparison with previous treatments 

The crucial difference between the treatment of the 
annealed disorder model given above and the one in Ref. 
4 is related to the occurrence of the coupling constant 
H( 2 K In perturbation theory, i.e., in an expansion in 
powers of q, the annealed disorder generates terms that 
have the structure of the last term on the right-hand 
side of Eq. (4.1b), except that they are not constrained 
to being diagonal in the replica index. There are two 
possible interpretations of such terms. (1) They could 
represent terms quadratic in Q that are not diagonal in 
replica space. This was the interpretation given in Ref. 
4. (2) They could present a new term linear in Q, which 



was not present in the original action. The term with 
coupling constant introduced in the present paper 
serves that purpose. By means of high-order perturba- 
tion theory one could in principle distinguish between 
these two possibilities, but this would be extremely cum- 
bersome. Let us instead argue on general structural and 
on physical grounds that the second interpretation is the 
correct one. 

First, we have argued in Sec. II that the annealed dis- 
order, since it gets averaged over at the level of the parti- 
tion function, should indeed be interpreted as an effective 
interaction between the electrons. As such, all involved 
degrees of freedom must occur with the same replica in- 
dex, and the generation of an interaction term (i.e., one 
quadratic in Q) for which this is not the case makes no 
physical sense. Terms quadratic in Q with more than 
one replica index are characteristic for quenched disor- 
der, and indeed the treatment of the annealed disorder 
in Ref. 4 was modeled after that of quenched magnetic 
disorder. As we have argued above, this is physically not 
plausible. 

Second, the appearance of a term with the structure 
of Aq \ Eq. (3.4), is plausible on physical grounds. The 
term with coupling constant in the NLaM repre- 

sents a frequency coupling with a microscopic time scale, 
on the order of an inverse Fermi energy (in units where 
h = 1). An interaction that is short-ranged in time does 
not add a new time scale to the problem. It therefore 
renormalizes but does not generate a new frequency 
coupling. An interaction that is long-ranged in time, 
on the other hand, does introduce a new time scale and 
hence a new frequency coupling. In the general case of 
a frequency dependent interaction with a continuum of 
time scales one would expect a frequency dependent cou- 
pling constant H whose scaling properties would have 
to be studied by means of a functional RG. In our sim- 
ple model where the annealed disorder is static, which 
means that the resulting effective interaction has an infi- 
nite range in time, one additional frequency coupling suf- 
fices, which is H( 2 \ The infinite time scale corresponds 
to a vanishing frequency scale, in accord with the discon- 
tinuous frequency dependence sgnfi in Eq. (3.4). In this 
context, we note that the term does not represent 
an inelastic lifetime. Rather, it is a true mass in the two- 
point propagators that is produced by the long-ranged 
in time interaction. This is analogous to the mass cor- 
responding to the plasmon pole that is produced by an 
interaction that is long-ranged in space. 

Third, the structure of the renormalization scheme 
used in Ref. 4 did not reflect the constraints discussed in 
Appendix B. This is only of minor concern if one neglects 
the coupling constant J^ 2,t > and uses one renormalization 
constant each for J s and Jt, as was done in Ref. 4. It be- 
comes crucial, however, in the presence of j( 2,t ), which 
forces the issue of how many renormalization constants 
are needed. 

Finally, the treatment of Ref. 4 led to results that 
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were not consistent with independent, very general, con- 
siderations. In particular, its prediction that dn/dfi is 
singularly renormalized, and critical at a MIT, contra- 
dicted one of the results of Rcf. 6. This point requires 
some explanation. The critical behavior predicted im- 
plies a nonanalytic dependence on the RG length scale, 
and hence a nonanalytic dependence on the wavenum- 
ber |q| in perturbation theory. In two-dimensions, this 
takes the form of a In |q| term in perturbation theory 
that is caused by the diffusive electron dynamics. In 
d > 2, these same integrals over diffusion poles lead to 
a |q| d ~ 2 dependence. 8 The predicted critical behavior of 
dn/dfi at the MIT, and the mechanism that causes it, 
therefore implies a nonanalytic wavenumber dependence 
of this susceptibility in the metallic phase. However, it 
was shown on general grounds in Ref. 6 that dn/dfi is an 
analytic function of the wavenumber for a large class of 
models, which includes the one under consideration here. 
This discrepancy prompted the current investigation, sec 
the discussion in Sec. I above. 



C. Conclusion, and Outlook 

In conclusion, we have found that the treatment in 
Ref. 4 of the electron problem in the presence of annealed 
disorder, in addition to quenched one, was not correct. 
The perturbation theory was correct, but the assump- 
tions made about the RG structure of the theory were 
not. This was the reason for the discrepancy between 
the explicit results found in Ref. 4 and later, more general 
considerations. 6 The current procedure, which considers 
the annealed disorder as an effective electron-electron in- 
teraction that is long-ranged in time, is physically and 
technically more convincing. It yields results that are 
consistent with all of the available information, and in 
particular with Ref. 6. Physically, the model of static, 
annealed magnetic disorder representing a type of local 
moments thus turns out to be less interesting than Ref. 
4 had given reason to believe. Instead of describing an 
unusual MIT, the model describes a variant of the ferro- 
magnetic transition of itinerant electrons that has been 
studied before. It is important to note that the same 
model with quenched instead of annealed magnetic dis- 
order is well known to contain a MIT in d = 2 + e. 23 This 
serves to underscore the fundamental physical difference 
between quenched and annealed disorder that we have 
stressed several times in this paper. 

We finally mention a possible consequence of our ob- 
servation, discussed in Sec. VB, that an electron-electron 
interaction with more than one time scale produces more 
than one frequency coupling in the NLctM, which in turn 
require additional renormalization constants. (In the 
present case, there was one additional time scale, infinity, 
one additional coupling, H^ 2 \ and one additional renor- 
malization constant.) At a MIT, the coupling constant 
that was denoted above by acquires a power-law 



frequency dependence. This is equivalent to saying that 
there are infinitely many time scales in the problem, and 
this raises doubts about the validity of renormalizing the 
action with just one renormalization constant for the fre- 
quency coupling. It is therefore possible that a complete 
description of the dynamics near a MIT would require a 
functional RG. A complete understanding of this prob- 
lem would also require a solution of the rcnormalizabil- 
ity problem for models of interacting electrons that is 
explained in Appendix B. 
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APPENDIX A: MOMENTUM-SHELL VERSUS 
FIELD-THEORETIC RENORMALIZATION 

In this appendix we motivate our choice of a field- 
theoretic formulation of the renormalization procedure. 

There exist two basic formulations of the RG, the field- 
theoretic one that originated in high-energy physics, 9 and 
Wilson's momentum-shell method, 17 which was invented 
for the study of critical points. After Wilson's break- 
through, it was shown that the field-theoretic method can 
also be applied to critical phenomena. 9 The relation be- 
tween these two formulations of the RG is complicated, 24 
but for our purposes we can restrict ourselves to a few 
basic features. 

In the Wilsonian method one renormalizes the Hamil- 
tonian or action itself, generating new interactions as one 
goes along, and checking all newly generated terms for 
their scale dimensions, and hence for their being RG rel- 
evant, irrelevant, or marginal. Irrelevant ones can be 
dropped, while relevant or marginal ones must be added 
to the model and included in a repetition of the renormal- 
ization process. In the field theoretic method, one renor- 
malizes specific propagators or vertex functions, and one 
needs to know from the outset how many renormalization 
constants are needed in order to make all of the vertex 
functions finite to all orders. 

For many models (e.g. for </> 4 -theory) there is only a 
small number of relevant or marginal terms. In these 
cases, the momentum-shell method is often preferred 
since it is physically more intuitive, and since it provides 
an explicit check for the generation of additional terms 
that must be kept. However, the NLctM does not belong 
to this class, as it has an infinite number of marginal 
terms in d = 2: In an expansion of Eq. (3.1) in powers 
of q, all terms are marginal. It is a priori unclear how 
the infinitely many coupling constants multiplying these 
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terms renormalize, although their bare values all coin- 
cide. The field-theoretic RG method proves that these 
coupling constants all renormalize the same way. 25 ' 9 This 
fixes the structure of the renormalized theory, and it then 
suffices to consider a small number of vertex functions 
in order to determine the renormalized theory explicitly. 
In the momentum-shell method, on the other hand, one 
needs to explicitly consider a large number of vertex func- 
tions or propagators (in principle infinitely many in the 
case of the NLcrM) in order to do the same. 

The same considerations apply to the terms in addi- 
tion to the NLcrM. Equations (3.7) add six coupling con- 
stants to the model. Within a momentum-shell RG, one 
would have to consider g 4 -vertices in order to determine 
how they renormalize. The field-theoretic method, on the 
other hand, allows us to argue that all of the J split into 
two pieces that pairwise renormalize in the same way, see 
Sec. IIIC and Appendix B. As a result, we need to ex- 
plicitly renormalize c/ 2 -vertices only. This is the reason 
why in this paper we choose the field-theoretic method 
over the momentum-shell one. 



APPENDIX B: INVARIANT DECOMPOSITION 
OF ANNEALED DISORDER TERMS 

In this appendix we recall the answer to the follow- 
ing question: Consider the NLcrM, Eq. (3.1), which is 
known to be renormalizablc in two-dimensions with two 
renormalization constants. 25 ' 9 . Now add to this action 
symmetry breaking operators. How does this affect the 
renormalizability, and how many additional renormaliza- 
tion constants are needed? 

For the case of operators that give some components 
of the basic field, Q(x) in our case, a mass (massive in- 
sertions), this question has been studied in detail. 26 If 
the NLcrM is invariant under transformations that form 
a symmetry group Q, then the operators in question must 
be expanded in a basis of irreducible representations of 
Q. All operators that belong to the same irreducible rep- 
resentation renormalize the same way, i.e., for each ir- 
reducible representation one additional renormalization 
constant is needed. 

In our case, it is most convenient to write the spin 
degrees of freedom explicitly, and consider the complex 

numbers Q"m,ij as tne matrix elements of Q. The 
NLcrM action is then invariant under unitary transfor- 
mations. We are interested in symmetry breaking oper- 
ators that are quadratic in Q. This case was first con- 
sidered by Pruiskcn. 18 There are two irreducible repre- 
sentations that correspond to symmetrized and antisym- 
metrized products of the Q. Any operator 

O = / "12,34 Qi2(x) Q 3 4(x) , (Bla) 

J 1234 

should thus be written as 



= 0++0^ , (Bib) 

with 

0± = ]r ( "12,34 [Ql2(x) Q 3 4(x) 

•* 12,34 

±Q 3 2(x)Qi 4 (x)] . (Blc) 

Here 1 = (ni,ai,ii), etc. + and require one renor- 
malization constant each, so two additional renormaliza- 
tion constants are needed to renormalize the NLcrM with 
arbitrary massive insertions of order Q 2 . 

A complication lies in the fact that in the present 
model, the coupling constants and multiply 

frequency dependent terms, and the frequency gets inte- 
grated over in perturbation theory. As a result, ratios 
of the J and H appear in perturbation theory, and the 
proof given in Rcfs. 26 docs not apply. This is true a 
fortiori in the case of an instantaneous electron-electron 
interaction, where the additional operators are not even 
massive insertions. Nevertheless, while no actual proof of 
renormalizability exists in this case, Ref. 19 has presented 
substantial evidence from perturbation theory that the 
model is still renormalizable with two additional renor- 
malization constants for the interaction. The same con- 
clusion is expected to hold in the annealed disorder case. 

APPENDIX C: COMPARISON WITH THE CASE 
OF AN INSTANTANEOUS INTERACTION 

In this appendix we compare the flow equations derived 
in Sec. IV B with those for the case of an instantaneous 
electron-electron interaction. 

In the instantaneous interaction case one has spin- 
singlet and spin-triplet interactions amplitudes K s and 
K t , that are analogous to j( 1,s ) and J^ 1 '*), respectively. 
The analog of j( 2,t ) does not exist. Instead of the two 
frequency couplings and there is only one cou- 
pling constant H, which is proportional to the specific 
heat coefficient, dn/d/j, and Xs are proportional to H+K s 
and H + K t , respectively. 1 

In the absence of magnetic impurities, a magnetic field, 
or spin-orbit scattering, K t flows towards large values, 
and after some transient behavior the one-loop flow equa- 
tions take the form 

^ = -eg + 0(g 3 ) , (Cla) 

^ = ^fct + 0(ff 2 ) , (Clb) 

^l = ^gk t + 0(g 2 ) , (Clc) 

§ = 5**?/* • (cm) 

A comparison with Eqs. (4.13) shows that the two be- 
haviors are very similar, except that in the instanta- 
neous interaction case k t flows to infinity much faster 
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than — k s . In particular, the conductivity and dn/d/i 
are not renormalized in either case (and neither is the 
density of states), while the magnetic susceptibility and 
the specific heat coefficient both diverge, albeit the lat- 
ter only logarithmically in the instantaneous interaction 
case. 10 Strictly at one-loop order, the physical interpre- 
tation of the RG flow was long considered not obvious, as 
has been stressed in the literature many times. However, 
the analysis given in Ref. 10, combined with the detailed 
discussion of the ferromagnetic transition in Ref. 20, has 
shown that the proper interpretation is in terms of a fer- 
romagnetic transition in d = 2 + e, as has been discussed 
in Sec. VA. 



APPENDIX D: PERTURBATION THEORY FOR 
THE FREE ENERGY 

Here we calculate the free energy in perturbation the- 
ory. This serves as a check on our scaling arguments for 
various observables in Sec. VA. 

To zeroth order in a loop expansion, the free energy 
density / is given by the saddle-point action. This 
yields free-electron values for all thermodynamic quanti- 
ties. The first correction, Af, is obtained by integrating 
over the fields q in Gaussian approximation. From Eqs. 
(4.1) we find 



A/ = A/ B + 3A/ t 



(Dla) 



Here 



A/ S , t = J s , t ^ d V (H^ + ^J Sjt) 

x- V / dun(w/T)P'' t (k ) w,T) , (Dlb) 
v ,. Jo 



with (cf. Eq. (4.2e)) 
P s>t (k,u;,T) 



k 2 -iGHMu + G(H<- 2 ) + t]J s , t )2nT ' 

(Die) 



The function 



n(x) = — coth( — ) 

'2 2' x 



(Did) 



serves as a convenient means for transforming the sum 
over Matsubara frequencies into a real-frequency integral, 
and we have used the familiar "charging formula" trick 
of integrating over the interaction constants in order to 
improve convergence. 
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